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This  paper  investigates  the  Proton  Exchange  Membrane  Fuel  Cell  (PEMFC)  Cold  Startup  problem  within 
the  specific  context  of  the  Plugin  Hybrid  Electric  Vehicles  (PHEV).  A  global  strategy  which  aims  at 
providing  an  efficient  method  to  minimize  the  energy  consumption  during  the  startup  of  a  PEMFC  is 
proposed.  The  overall  control  system  is  based  on  a  supervisory  architecture  in  which  the  Energy 
Management  System  (EMS)  plays  the  role  of  the  power  flow  supervisor.  The  EMS  estimates  in  advance, 
the  time  to  start  the  fuel  cell  (FC)  based  upon  the  battery  energy  usage  during  the  trip.  Given  this 
estimation  and  the  amount  of  additional  energy  required,  the  fuel  cell  temperature  management  strategy 
computes  the  most  appropriate  time  to  start  heating  the  stack  in  order  to  reduce  heat  loss  through  the 
natural  convection.  As  the  cell  temperature  rises,  the  PEMFC  is  started  and  the  reaction  heat  is  used  as 
a  self-heating  power  source  to  further  increase  the  stack  temperature.  A  time  optimal  self-heating 
approach  based  on  the  Pontryagin  minimum  principle  is  proposed  and  tested.  The  experimental 
results  have  shown  that  the  proposed  approach  is  efficient  and  can  be  implemented  in  real-time  on  FC- 
PHEVs. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  hydrogen  fuel  cell  is  emerging  as  one  of  the  best  power 
sources  for  the  sustainable  transportation  due  to  its  low  green¬ 
house  gas  emission  and  its  high  power  density  [1,2].  Different  types 
of  fuel  cell  technologies  are  being  investigated  for  Hybrid  Electric 
Vehicle  (HEV).  However,  the  Proton  Exchange  Membrane  Fuel  Cell 
(PEMFC)  is  one  of  the  most  selected  power  sources  for  electric 
mobility  [3].  In  some  vehicle  applications  involving  fuel  cells,  the 
batteries  are  being  used  primarily  to  capture  the  regenerative 
power,  secondary  to  absorb  the  rapid  and  high  power  transient 
originating  from  the  vehicle  dynamics  and  finally  to  provide  power 
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to  the  vehicle  motion  [4].  These  vehicles  are  designated  as  Fuel  Cell 
Plug-in  Hybrid  Electric  Vehicle  (FC-PHEV).  The  energy  efficiency  of 
an  FC-PHEV  is  thoroughly  related  to  the  power  flow  management 
between  all  components  that  contribute  to  the  vehicle  propulsion. 
So,  a  lot  of  work  is  being  done  on  the  specific  problem  of  the  on¬ 
board  energy  management  and  power  splitting  [5]. 

When  the  fuel  cell  is  operating,  there  is  a  production  of  water 
and  water  vapor.  Although  a  systematic  purge  is  performed  during 
the  PEMFC  shutdown  process,  the  remaining  water  vapor  could 
frost  when  the  temperature  drops  below  0  °C  (in  regions  with  long 
and  very  cold  winter  seasons).  To  preserve  the  fuel  cell  life, 
a  specific  attention  is  required  before  restarting  the  fuel  cell  during 
the  periods  of  sub-freezing  and  freezing  ambient  temperatures 
[6-8].  Thus,  Ref.  [9]  reported  that  an  irreversible  degradation  in  the 
fuel  cell  components  occurs  during  sub-freezing  and  freezing 
operating  temperatures.  The  volume  expansion  of  the  cell  materials 
in  the  freezing  conditions  is  one  of  the  reasons  for  this  degradation 
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[10].  Among  the  proposed  recommendations  in  Refs.  [9,10]  is  the 
water  purge  operation  during  the  fuel  cell  shutdown.  The  research 
reported  in  Refs.  [11,12]  showed  that  heating  the  stack  has  a  great 
positive  impact  on  the  fuel  cell  transient  thermal  model.  In  Ref.  [13], 
the  catalytic  oxygen/hydrogen  reaction  heat  is  used  to  increase  the 
stack  temperature. 

In  physics  and  in  electrochemistry,  several  studies  regarding  the 
stack  self-heating  when  the  stack  temperature  belongs  to  a  sub¬ 
freezing  range,  have  been  reported.  The  cell  membrane  water 
content  and  the  ice  formation  in  the  cathode  catalyst  layer  are 
extensively  studied  and  identified  as  the  main  performance 
limiting  factors.  In  Ref.  [21  ],  an  analytical  model  describing  the  heat 
balance,  the  ice  formation  in  the  catalyst  layer,  and  the  water 
transport  characteristic  throughout  a  fuel  cell  at  very  low  temper¬ 
atures  is  presented.  In  Ref.  [22],  numerical  studies  were  conducted 
using  a  monotonically  increasing  current  for  a  fuel  cell  cold  start 
from  -30  °C.  It  was  shown  that  using  a  moderate  current  increasing 
rate  or  using  a  large  initial  current  density  in  combination  with 
a  small  current  increasing  rate,  can  lead  to  a  successful  stack  self¬ 
start.  In  the  experimental  study  presented  in  Ref.  [23],  the  ice 
accumulation  in  a  frozen  PEMFC  was  measured  at  sub-freezing 
temperatures  and  the  time  to  voltage  failure  for  various  current 
densities  under  isothermal  operation  at  -20  °C  was  also  measured. 

Although  these  reported  studies  demonstrated  that  it  is  possible 
to  achieve  a  fuel  cell  self-start  from  a  sub-freezing  temperature,  we 
propose  in  this  paper  a  different  approach  which  focuses  only  on 
the  stack  thermodynamic  and  electric  aspects.  This  approach 
allows  us  to  consider  the  PEMFC  global  electrothermal  behavior 
without  focusing  on  the  electrochemistry  detail  equations.  In 
addition,  the  electrothermal  behavior  model  permits  a  simple  and 
effective  formulation  of  the  global  optimization  problem  that  can 
be  solved  in  real-time  in  the  context  of  the  FC-PHEV.  To  avoid 
complicated  water  management  issues  and  potential  irreversible 
damages  in  the  fuel  cell  components,  an  external  heating  system  is 
used  to  heat  the  stack  and  let  the  cell  temperature  to  belong  to  its 
normal  operating  range  from  a  sub-freezing  starting  point. 

Two  different  methods  have  been  reported  for  the  fuel  cell 
heating  [7].  The  internal  heating  which  is  the  first  method,  consists 
in  providing  an  electric  current  to  a  resistance  embedded  directly 
into  the  stack.  The  heat  generated  from  the  resistance  Joule  effect  is 
transferred  to  the  stack  structure.  Often,  this  method  is  used  by  the 
fuel  cell  manufacturers.  The  second  method  is  related  to  the 
external  fluid  based  warm-up  system  which  uses  a  heat  exchanger 
and  an  antifreeze  fluid  (mixture  of  glycol  and  water).  The  exchanger 
allows  the  heat  transfer  from  the  fluid  to  the  stack  components. 
This  method  is  more  accessible  to  our  research  than  the  previous 
one.  However,  the  main  problem  is  to  determine  the  heating  energy 
(heating  power  and  duration)  that  will  allow  rising  the  stack 
temperature  whilst  minimizing  the  heat  loss  through  the  natural 
convection  with  the  environment. 

This  paper  investigates  the  PEMFC  temperature  management 
for  the  vehicle  application,  especially  for  the  FC-PHEV  vehicles.  On 
such  a  vehicle,  the  batteries  have  the  energy  to  propel  the  vehicle 
over  a  portion  of  the  travel  distance  and  the  fuel  cell,  which  plays 
the  role  of  range  extender,  is  started  to  recharge  the  batteries 
through  a  serial  topology  [14]. 

We  assume  that  the  fuel  cell  has  an  appropriate  heat  exchanger 
embedded  into  the  cells  to  allow  the  warm-up  process  using  an 
external  fluid  heating  system.  Two  contributions  are  presented:  the 
first  one  is  related  to  the  design  and  to  the  experimental  validation 
of  a  two  layer  supervisory  architecture  for  FC-PHEV  energy 
management.  The  second  contribution  consists  in  designing 
a  global  strategy  for  the  fuel  cold  startup  by  considering  the  cells 
initial  temperature,  the  ambient  temperature,  the  thermal  oper¬ 
ating  range  of  the  cells  and  the  amount  of  energy  available  for 


heating  the  stack.  The  vehicle  Energy  Management  System  (EMS)  is 
part  of  the  upper  layer  whereas  the  global  strategy  for  the  fuel  cell 
cold  startup  belongs  to  the  lower  layer.  Knowing  the  trip  length  and 
the  average  power  requirement  for  the  vehicle  motion,  this  EMS 
estimates  the  amount  of  fuel  cell  energy  required  to  complete  the 
trip.  It  also  estimates  in  advance  the  appropriate  time  to  start  the 
fuel  cell.  The  global  strategy  for  the  fuel  cell  cold  startup  uses  these 
two  data  as  constraints  and  provides  the  appropriate  heating 
timing  that  will  allow  the  fuel  cell  to  comply  with  them. 

The  rest  of  the  paper  is  organized  in  7  sections.  The  FC-PHEV 
supervisory  architecture  is  presented  in  section  2.  The  cold  start 
problem  overview  and  the  startup  strategy  are  discussed  in 
sections  3  and  4,  respectively.  The  self-heating  method  based  on 
the  Pontryagin  minimum  principle  is  investigated  in  section  5, 
whereas  the  high  temperature  regulation  and  the  experimental 
results  are  discussed  in  sections  6  and  7,  respectively.  Finally  the 
conclusion  is  presented  in  section  8. 


2.  FC-PHEV  supervisory  architecture 

Fig.  1  represents  a  typical  FC-PHEV  serial  topology  where  the 
solid  and  dashed  arrows  represent  the  power  flow  and  the  control 
signal,  respectively.  The  driver  power  demand  is  sent  to  the  Energy 
Management  System  (EMS)  which  generated  the  corresponding 
power  command  Pm  for  Electric  Propulsion  System.  The  bi¬ 
directional  converter  allows  the  Battery  Pack  to  be  charged  and 
discharged  and  Pb  represents  the  power  flowing  through  this 
converter.  In  this  architecture  the  EMS  supervises  the  Fuel  Cell 
Power  Controller  and  the  battery  pack  energy  E b.  The  power  from 
the  Fuel  Cell  PfC  is  used  to  recharge  the  batteries. 

To  prevent  the  Battery  Pack  to  be  over  discharged,  a  minimum 
energy  threshold  E^m  is  usually  defined  [15].  So  the  Battery  Pack  is 
operating  in  a  depleting  mode  until  Eb  reaches  E™in.  Knowing  the 
initial  Battery  Pack  energy  E{],  the  amount  of  energy  available  for 
the  trip  and  the  stack  heating  is  represented  by: 

El  =  £0-E™n-£ma*  (1) 

where  E™x  is  the  maximum  energy  that  can  be  used  to  heat  the 
stack. 

We  assume  that  the  whole  trip  duration  tt d  and  the  average 
traction  power  Pm  are  known;  thus,  it  is  possible  to  approximate 
the  total  energy  required  to  drive  the  vehicle  to  the  destination 
using  the  following  equation: 
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Fig.  1.  Fuel  cell  hybrid  electric  Vehicle:  the  arrows  with  solid  lines  represent  energy 
flows;  the  arrows  with  dashed  lines  represent  control  and  measured  signals. 
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Etr  —  Emttd  (2) 

If  (Eg  <  Et r),  the  additional  energy  required  for  the  Fuel  Cell  is 
therefore  estimated  using  the  following  equation: 

Efc  =  Etr-Eg  (3) 

As  the  battery  energy  Eb  reaches  E^in,  the  fuel  cell  is  started  to 
recharge  batteries.  Since  the  fuel  cell  needs  to  be  warm-up  before 
being  started,  the  EMS  must  anticipate  when  the  batteries 
minimum  energy  E^m  will  be  attained.  If  Pm  is  perfectly  predict¬ 
able,  then  the  time  tst  to  reach  E^in  is  obtained  by  solving  the 
following  equation: 

tst 

E°b  -  J  Pm(t)dt  -  E£™  =  £™n  (4) 

0 

where  Eg  is  the  Battery  pack  initial  energy  and  where  tst  >  0. 

In  practice,  the  perfect  knowledge  of  Pm  which  depends  on  the 
value  of  Pu  is  difficult  to  obtain  because  the  vehicle  driver  command 
is  related  to  the  unpredictable  factors  such  as  the  traffic,  the  road 
conditions,  etc.  So,  an  approximation  tst  of  tst  is  required. 


environment.  The  operation  of  this  system  will  help  the  stack 
temperature  to  increase  from  its  sub-freezing  point  to  the  Tf™in. 

As  the  stack  temperature  TfC  reaches  Tf™m,  the  fuel  cell  can  be 
started  safely.  The  electrochemical  model  of  the  PEMFC  shows  that 
the  efficiency  is  related  to  several  parameters  such  as  the  temper¬ 
ature,  the  air  humidity  ratio  and  the  air  stoichiometry.  Since  we  are 
mainly  interested  in  the  temperature  effect  on  the  PEMFC  work 
point,  several  studies  and  analyses  have  demonstrated  that  by 
increasing  the  temperature,  the  fuel  cell  efficiency  is  improved  for 
high  power  production  [24,25]. 

When  TfC  reaches  Tf™in,  it  is  not  efficient  to  use  the  battery 
energy  to  heat  the  fuel  cells  in  order  to  raise  its  temperature 
quickly.  Instead,  the  heat  from  the  exothermic  reaction  in  the  fuel 
cell  is  used  to  increase  the  stack  temperature.  The  second  problem 
is  to  find  an  optimal  fuel  cell  power  production  profile  in  order  to 
increase  as  quickly  as  possible  the  stack  temperature  without 
exceeding  the  prescribed  maximum  temperature  Tf™ax. 

When  the  stack  temperature  reaches  Tg}ax,  a  cooling  process 
must  be  started  to  prevent  overheating  cells.  The  PEMFC  high 
temperature  regulation  has  been  solved  in  the  literature  and 
numerous  methods  have  been  proposed.  Therefore,  it  will  not  be 
analyzed  in  this  paper. 


2  A.  Estimation  of  tst 

Assume  that  at  time  t  <  ts t,  the  past  values  of  Pm(t)  are  known.  So 
the  time  to  reach  Eg1111  is  given  by: 

t  tst 

El-  j  Pm(T)dT  -  J  Pm(T)dT  -  E£ax  =  E™n  (5) 

0  t 

Since  Pm  is  the  average  power  obtained  with  the  past  values  of  Pm, 
Eq.  (5)  becomes: 

t 

El-  j  Pm(t)dT  -  Pm  (tst  -  t)  -  £™x  =  EJT  (6) 

0 

Elence,  ts t  is  obtained  using  the  following  equation: 

t 

El~  J  Pm(t)dT-E™n-E^x 

fst  =  t  + - 5 - ~ -  (7) 

Pm 

3.  Cold  start  problem  overview 


4.  Cold  startup  strategy 

4 A.  Stack  warm-up  strategy 

In  this  section  we  propose  a  method  that  will  provide  the  most 
appropriate  time  tst,he  to  start  the  warm-up  system  to  increase  the 
fuel  cell  stack  temperature  using  the  following  assumptions: 

•  the  time  to  start  the  PEMFC  fst  is  given  by  the  EMS; 

•  the  amount  of  energy  to  produce  is  also  given; 

•  the  heating  fluid  is  incompressible; 

•  the  PEMFC  model  is  represented  by  the  heating  thermal  mass 
mfc  and  the  average  specific  heat  capacity  Cfc; 

•  he  heating  power  of  the  warm-up  system  is  fixed  and  known 
beforehand. 

Fig.  2  represents  the  warm-up  system  connected  to  the  PEMFC 
stack. 

The  first  order  heat  transfer  model  of  the  system  is  represented 
by  the  following  system  of  differential  linear  equations: 

tnwCyvEw  —  Phe  Q-vvhe  (g) 

tftfcCfcPfc  =  —  Q-amb  +  Q.whe 


The  target  of  the  cold  startup  strategy  is  related  to  the  Flybrid 
Electrical  Vehicles  (HEVs)  with  a  PEMFC  as  a  range  extender  power 
source.  We  assume  that  the  vehicle  EMS  provides  the  estimated 
startup  time  tst  and  the  estimated  fuel  cell  energy  production  Efc- 
Furthermore,  we  assume  that  the  initial  temperature  Tfc  of  the 
PEMFC  as  well  as  the  ambient  temperature  Tam b  are  also  known 
beforehand.  If  and  Tamb  are  in  the  sub-freezing  temperature 
range,  it  is  recommended  to  warm-up  the  fuel  cell  above  the 
freezing  temperature  in  order  to  avoid  complex  fluid  management 
issues,  ice  build-up  and  to  preserve  the  cell’s  life  [9,10]. 

Indeed,  the  PEMFC  have  an  operating  temperature  range  (Tf™in  < 
Tfc  <  Tf™ax)  which  depends  on  the  type  of  technology  and  the  factory 
specifications.  Since  the  amount  of  energy  required  to  warm-up  the 
fuel  cell  is  limited  to  E[^ax,  the  first  problem  to  solve  is  to  determine 
the  appropriate  time  to  start  the  fuel  cell  warm-up  system  in  order 
to  minimize  the  heat  transfer  from  the  fuel  cell  structure  to  the 
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Fig.  2.  Warm-up  system  and  PEMFC. 
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where  mw,  mfC,  cw,  CfD  Tw,  Phe,  Qamp,  and  QWhe  represent  the  heating 
fluid  mass,  the  stack  mass,  the  average  specific  heat  coefficient  of 
the  heating  fluid,  the  average  specific  heat  coefficient  of  the  fuel 
cell,  the  heating  fluid  average  temperature,  the  heating  element 
power,  the  stack-environment  heat  transfer  rate  (see  Eq.  (9))  and 
the  heating  element  heat  transfer  rate  (see  Eq.  (10)),  respectively. 

Q.amb  —  hnc$nc  —  ^ambj  (9) 

where  hno  Sno  Tamp  and  TfC  represent  the  natural  convective  heat 
transfer  coefficient,  the  effective  heat  transfer  surface  (between  the 
stack  and  the  environment),  the  ambient  temperature  and  the  fuel 
cell  temperature,  respectively. 

In  the  considered  open  thermodynamic  system,  the  heat 
transferred  by  the  heating  fluid  is  represented  by: 

Q-whe  =  thwCvv  (Tw.out  —  fw,in  )  (10) 


where  mw,  Tw> jn  and  Tw>out  represent  the  heating  fluid  mass  flow 
rate,  the  temperatures  of  the  heating  fluid  that  is  entering  and 
exiting  the  stack,  respectively. 

The  expression  for  the  outlet  fluid  temperature  is  given  by  [20]: 


=  Tfr 


(lw,in  Tfc)  ( 


frcw^cw 

rhwcw 


(11) 


where  hcw  and  5CW  represent  the  forced  convection  heat  transfer 
coefficient  and  the  heat  transfer  surface,  respectively. 

Substituting  TWi0Ut  in  Eq.  (10)  by  its  expression  (11)  yields: 


Q-whe  —  kw  ^lw,in 

where: 


(12) 


kw  —  mwcw  ^1  —  e 


frcwScw 

rhwcw 


The  solution  of  the  system  for  differential  Eq.  (8)  is  given  by: 


7>c(f)  =  Tf”  +  (Tf“  -  Tf“  -  c)e-r'r  +  Ce~r2t 


(13) 


where: 


rf  oo  _  7-  j  *■  he 

J  fc  —  J  amb  “I"  u  <- 
u  nc^n< 


(14) 


and  where: 


C  = 


/<wTo  ,  kw  To  ,  hncSncT  tq  r  T  oo  ^ 

1  w  "I"  7^“  J  fc  "f  7^  J  amb  —  1  fc'  2  —  1  fc  '  1 

'-fc  '-W  '-fc 

(n  -  r2) 


(15) 


and  where: 


1  /  hnc$nc 
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Qc  Qc  Cw 


1 
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(16) 
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(17) 


And  where  Cfc  =  rrifCCfC  and  Cw  =  mwcw  are  the  heat  capacity  of  the 
PEMFC  and  of  the  heating  fluid  respectively. 

Consider  that  the  ambient  temperature  is  lower  than  Tg}m  (the 
minimum  prescribed  temperature  for  the  fuel  cell  operation).  The 
role  of  the  warm-up  system  is  to  increase  TfC  up  to  Tg}m  by  heating 
the  stack.  However,  if  the  system  starts  heating  too  earlier,  a  part  of 
the  energy  is  lost  into  the  environment  through  the  natural 
convection.  So  in  the  next  section,  we  propose  a  simple  method 
that  will  allow  determining  the  most  appropriate  time  to  start 
heating  the  stack. 

4.2.  Time  to  start  heating  the  fuel  cell  stack 

To  minimize  the  heat  transfer  between  the  stack  and  the  envi¬ 
ronment  during  the  warm-up  process,  the  heating  duration  must 
be  short  and  adequate  to  bring  the  fuel  cell  temperature  close  the 
T™in.  So,  given  the  time  to  start  the  fuel  cell  ts t»  the  shortest  heating 
duration  the  is  computed  using  Eq.  (13).  In  practice,  we  have  found 
that  the  following  approximation  can  adequately  represent  the 
heating  duration: 


the 


Cfcw 

hnc^nc 


In 


knc^nc  y 

jmm  _  j 

1  fo  1  amb 

)  —  Qie 

hnc$nc 

(Tl  -  Qmb) 

—  Phe  - 

(18) 


where: 


Cfow  =  mwCw  +  mfcCfo  (19) 

So,  by  taking  into  account  the  estimated  time  to  start  the  fuel  cell 
?st,  we  obtained  the  expression  of  the  time  to  start  heating  the  stack 
tst,he  as  followed: 

tst,he  ~  tst  —  the  (20) 

The  total  energy  used  from  the  batteries  by  the  warm-up  system 
is  therefore  approximated  as  follows: 

tst 

t^he  =  J  ^he^t  =  Phe(^tst  ~  ^st,he)  (21) 

fst,he 


5.  Stack  self  heating 

As  the  stack  temperature  reaches  the  minimum  prescribed 
temperature  T™in,  the  PEMFC  can  be  started.  The  efficiency  which 
depends  on  the  cell  temperature  is  low.  Increasing  the  temperature 
by  providing  external  heat  can  help  to  improve  this  efficiency  and 
avoid  ice  formation  near  the  freezing  temperature.  However,  this 
method  is  inefficient  since  the  electric  energy  used  to  accelerate  the 
stack  warm-up  is  no  longer  available  for  the  vehicle  traction.  Since 
the  electricity  production  in  the  PEMFC  is  exothermic,  the  heat 
from  the  reaction  can  be  used  as  a  self-heating  energy  to  increase 
the  stack  temperature.  In  this  section,  we  propose  a  method  to 
determine  the  optimal  level  of  the  fuel  cell  power  production  in 
order  to  increase  quickly  the  stack  temperature  and  avoid  ice 
formation. 

5.1.  Optimal  current  profile 

5.1.1.  Problem  formulation 

As  the  fuel  cell  temperature  reaches  Tf™in,  the  PEMFC  is  started 
and  the  warm-up  system  is  stopped.  The  evolution  of  the  PEMFC 
temperature  depends  on  the  amount  of  heat  provided  by  the  two 
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heat  sources:  the  mass  species  heat  rate  Q.mst  and  the  reaction  heat 
rate  QYh.  The  heat  introduced  by  the  reaction  depends  on  the 
PEMFC  efficiency,  and  it  is  affected  by  several  operating  parameters 
such  as  the  load  current,  the  cell  temperature,  the  membrane  water 
content,  air  stoichiometry  ratio,  the  anode  and  the  cathode  pres¬ 
sures,  etc.  The  dynamics  of  the  fuel  cell  average  temperature  are 
represented  by  the  following  first-order  differential  equation: 

mfcCfc^  =  Qa  +  Qrh  -  Qamb  (22) 

where  Q.amb  is  represented  by  Eq.  (9);  Q.rh  is  the  reaction  heat 
transfer  rate  (see  Eq.  (25));  Qa  is  the  air  flow  heat  transfer  rate  (see 
Eq.  (23)).  The  hydrogen  heat  transfer  rate  is  neglected  in  this  model. 

As  in  the  case  of  forced  convection  (see  Eqs.  (10)— (12)),  the  heat 
transfer  rate  between  the  air  and  the  stack  is: 

Qa  «  /<a(Tain-Tfc)  (23) 

where: 


So  we  formulate  a  minimum  time  optimization  problem  as  follows: 
find  the  optimal  heat  rate  u(t)  that  minimizes  the  cost  function 
represented  by  Eq.  (31)  under  the  constraint  of  the  stack  temper¬ 
ature  dynamics  represented  by  Eq.  (27). 

£fc,min 

J=  j  dt  (31) 

tst 

where  tfC,min  represents  the  minimum  time  it  takes  for  the  fuel  cell 
temperature  to  reach  the  maximum  prescribed  temperature  T^ax. 
1st  is  obtained  by  Eq.  (7). 

5.1.2.  Solving  the  time  optimal  problem 

The  Hamiltonian  of  the  system  is  represented  by  [16]: 

jr(Tfc,  w)  =  1 +A(t)  (A7fc(r)  +Bu(t))  (32) 

where  A  represents  the  system  costate  variable  which  differential 
equation  is  given  by: 


/<a 


ha$a\ 

rnaca) 


(24) 


0jr(7fC,A,u,t) 

Wc 


(33) 


And  where  ma,  ca,  ha,  Sa  and  Ta  in  represent  the  air  flow  rate,  the 
specific  heat  coefficient  of  the  air,  the  forced  convection  heat 
transfer  coefficient  of  the  air,  the  heat  transfer  surface  between  the 
air  and  the  stack  and  the  inlet  air  temperature,  respectively.  The 
reaction  heat  rate  is  represented  by  [17]: 


Differentiating  Eq.  (32)  and  using  Eq.  (33)  yields: 

-A(f)  =  Ak(t)  (34) 

Hence,  the  solution  of  this  differential  equation  is  given  by: 


Qrh 


(  AHrh 

v  2  F 


-  v, 


cell  )ifcJ 


IfrN 


(25)  A(t)  =  A  (tfcmin)  e_A(tfcmin_t) 


(35) 


where  AHrh,  F,  Vceii,  kc  and  N  represent  the  enthalpy  of  the  forma¬ 
tion  [18],  the  Faraday  constant,  the  cell  voltage  as  represented  in 
Ref.  [19],  the  stack  output  current  and  the  number  of  cells. 

The  cell  voltage  is  given  by: 


Since  tfC,min  is  free,  dtfC>min  ^  0.  The  transversality  condition 
becomes: 

jr(Tfc,A,u,tfcmin)  =  0  (36) 


K:ell  —  Vocv  —  K)hm  —  ^act  —  ^con 


(26)  This  equation  is  equivalent  to: 


where  \70CV,  V0hm,  ^act  and  VCon  represent  the  open-circuit  voltage  or 
Nernst  voltage,  the  ohmic  loss,  the  activation  loss  and  the 
concentration  loss. 

Linearizing  Eq.  (22)  yields: 


1  + 

,min  )  (^fc  (tfc  ,min  )  +Bu(tfc 

,min  ))  =  0  (37) 

Since  we  would  like  to  have  Tfc(tfc  min)  =  T^ax,  the  costate  value 
A(tfC,min)  is  computed  as  follows: 


Tfc  =  ATfc  +  Bu  (t) 


where: 

^  _  ka  +  hCn$cn 
mfcCfc 

and  where: 

B  =  —l— 

mfccfc 

and  where: 


(27) 


(28) 


(29) 


U(f)  ~  Q.rh  +  ^cn5CnFamb  +  kaTa  in  (30) 

Given  the  temperature  dynamics  represented  by  Eq.  (27),  the 
objective  is  to  find  the  current  profile  IfC  (see  Eqs.  (25)  and  (30))  that 
will  quickly  increase  the  stack  temperature  from  Tj?}m  to  Tf^ax.  This 
current  must  not  exceed  the  maximum  prescribed  current  /™ax  >  0. 


^(ffc,min)  —  /  \  (38) 

V  )  AT™ax  +  Bu  (tfc  min) 

If  a  minimum  time  optimal  control  u*  exists,  then  the  Pontryagin 
minimum  principle  states  that  [16]: 

jf(r;c,AV,t)  <jr(r;c,x*,u,t)  (39) 

where  T^,  X*  and  u*  represent  the  optimal  fuel  cell  temperature 
profile,  the  optimal  costate  profile  and  the  optimal  heat  rate, 
respectively. 

Using  the  Hamiltonian  Eq.  (32),  the  expression  (39)  becomes: 
A*  (fj Bu *  (7)  <  A*  (fj  Bu  (fj  (40) 

Hence,  the  optimal  heat  rate  u*(t)  is  obtained  by  the  following 
equation  [16]: 

u*(t)  s=  -sgn (BA (t))umax  (41) 

where: 
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f  1,  if  BA(t)>0 

sgn (BX(t))  =  <  indeterminate,  if  BX(t)  =  0  (42) 

[-1,  if  BX(t)  <  0 

and  where  umax  is  the  maximum  heat  rate  that  can  be  applied  to  the 
stack. 

Eq.  (42)  shows  that  the  optimal  heat  rate  will  always  take  the 
extremum  values  of  u.  Hence  u*{t )  will  take  only  two  values: 


leads  to  more  heat  generation  than  electric  production  (startup 
stage).  When  the  temperature  reaches  the  maximum  prescribed, 
the  efficiency  is  high  which  reduces  the  proportion  of  the  heat 
versus  the  electricity. 

7.  Experimental  validation  and  discussions 

7.1.  Experimental  goal 


•  u*{t)  =  umax  if  (BA(t))  is  negative; 

•  u*{t)  =  -umax  if  (BA(t))  is  positive. 

In  the  case  of  A(t)  =  0,  there  is  a  singularity.  Fortunately,  this  case 
will  not  occur  if  umax  is  very  high  as  demonstrated  below.  Assume 
that  u*(tfc,min)  =  umax;  the  final  value  of  the  optimal  costate  is  given 
by: 


The  experiments  aim  at  validating  the  following  methods: 

the  estimation  of  the  Energy  Management  System  of  the  time 
to  start  the  fuel  cell; 

the  estimation  of  the  time  to  start  heating  the  stack; 

the  fuel  cell  thermal  model  during  the  stack  self-heating  stage; 

the  time  minimal  solution  approach  proposed  in  this  paper. 


A 


1 

AT™  +  Bumax 


(43) 


So,  the  costate  profile  which  is  derived  from  Eq.  (35)  becomes: 


X(t)  = 


1 

AT™*  +  Buma!<e 


^(tf  c,min  t) 


(44) 


If  there  is  enough  heating  power  to  warm-up  the  stack  then: 


jjmax 


Tmax 

lfc 


(45) 


In  this  case,  A(t)  <  0  which  confirms  that  i/(tfC,min)  =  umax.  In 
addition,  A(t)  cannot  be  equal  to  0.  Clearly,  there  is  no  singularity 
with  this  control  law  under  the  condition  (45).  The  minimum  time 
optimal  solution  is  therefore  to  apply  umax.  The  maximum  reaction 
heat  rate  Qrh  which  corresponds  to  umax  can  be  computed  using 
Eq.  (30).  Then,  the  optimal  current  profile  is  obtained  by  Eq.  (25). 


6.  High  temperature  regulation 

When  the  fuel  cell  temperature  reaches  the  maximum  prescribed 
temperature  r^ax,  it  should  be  regulated  in  order  to  avoid  over¬ 
heating  the  stack.  The  high  temperature  regulation  has  been  exten¬ 
sively  studied  in  the  literature.  In  this  paper,  a  simple  PID  regulator 
was  used  to  maintain  the  fuel  cell  temperature  close  to  Tg13*. 

Fig.  3  illustrates  the  different  energies  involved  in  the  PEMFC 
during  the  course  of  a  driving  cycle.  During  the  first  stage  (pre¬ 
startup  or  warm-up  phase),  external  heating  power  is  provided  to 
the  stack.  During  this  pre-startup  phase  there  is  no  power 
production  from  the  PEMFC.  As  the  stack  temperature  is  rising,  this 
external  power  is  removed  and  the  fuel  cell  is  started  to  begin 
producing  heat  (blue  color  zone)  and  electric  power  (shaded  green 
zone).  However,  the  low  efficiency  at  the  beginning  of  its  operation 


Total  energy  distribution  in  a  driving  cycle 


Pre  -  startup  startup  optimal  operating  temperature 


7.2.  Experimental  setup 

A  specific  test  bench  was  built  with  a  small  commercial  500  W 
PEMFC  which  picture  is  shown  in  Fig.  4.  The  experimental  param¬ 
eters  are  presented  in  Table  1.  The  complete  setup  is  represented  in 
Fig.  5.  In  this  figure,  we  have  an  air  supply  circuit,  a  hydrogen  supply 
circuit  operating  in  dead-end  conditions  and  a  fluid  cooling  circuit. 

The  air  supply  circuit  includes  a  compressor,  a  mass  flow  meter, 
an  air  pressure  sensor,  and  a  solenoid  valve.  To  control  the  air 
relative  humidity  and  the  inlet  air  temperature,  we  use  a  humidifier 
and  a  thermostatic  bath.  The  hydrogen  circuit  includes  a  hydrogen 
tank  a  solenoid  which  is  periodically  opened  to  purge  water  on  the 
cathode  side  of  the  fuel  cell  stack. 

The  fluid  cooling  circuit  includes  a  thermostatic  bath  (see  Fig.  6) 
which  is  controlled  using  a  personal  computer.  An  electronic  vari¬ 
able  load,  controlled  by  the  computer,  is  connected  to  the  output  of 


Fig.  3.  Energy  distribution  during  a  driving  cycle. 


Fig.  4.  Hyteon  500  W  PEMFC  system. 


Table  1 

Experimental  parameters. 
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Parameters 

Values 

Stack 

N 

23  (Cell) 

rrifc 

12  kg 

Cfc 

770  (J  kg"1  K_1) 

Natural  convection 

hnc 

5  (W  K1  irr2) 

$nc 

0.24  (m2) 

Air  forced  convention 

Ca 

1.0  (kj  kg"1  K"1) 

ha 

90  (W  K1  m-2) 

5a 

0.4968  (m2) 

ma 

3.57  x  10-7Aa/fCN(kg  s"1) 

Aa 

2.4 

Water  forced  convection 

cw 

4200  (kj  kg"1  K-1) 

mw 

2  (kg) 

hew 

34  (W  K1  m-2) 

Scw 

0.5962  (m2) 

rhw 

0.055  (kgs"1) 

the  fuel  cell  stack.  This  variable  load  is  used  to  emulate  the  FC-PHEV 
power  demand  during  the  experiments. 

The  computer  hosts  all  the  Labview  programs  relative  to  the 
methods  described  throughout  this  paper.  The  main  user  interface 
is  represented  in  Fig.  7.  The  hardware  for  the  signal  acquisition  and 
the  system  control  is  built  with  the  National  Instrument  devices. 

Our  setup  was  unable  to  run  the  experiments  with  low 
temperature  conditions.  So  we  decide  to  perform  the  tests  with 
a  temperature  range  of  [293.15  K,  323.15  K].  This  range  is  sufficient 
to  test  the  proposed  methods  and  to  validate  the  overall  PEMFC 
temperature  management  strategy. 

7.3.  Test  scenario 

7 3 A.  Driving  cycle  and  FC-PHEV  power  demand 

A  driving  cycle  obtained  when  an  experimental  electric  vehicle 
was  moving  on  the  campus  of  “Universite  du  Quebec  a  Trois-Riv- 
ieres”  is  used  for  the  test  and  the  validation  of  the  proposed 
methods.  Knowing  the  vehicle  parameters  and  the  driving  cycle, 
we  computed  the  power  demand  Pm.  A  down  scale  coefficient  is 
applied  to  the  calculated  Pm  in  order  to  match  the  fuel  cell  power 
production  which  maximum  value  is  set  to  500  W.  This  power 
demand  is  represented  in  Fig.  8. 

7.3.2.  Energy  management  system  (EMS)  parameters 

The  battery  energy  is  represented  by  the  State-Of-Charge  (SOC) 
which  is  the  ratio  between  Eb  and  Ef].  At  the  beginning  of  the 
experiment,  the  State-Of-Charge  (SOC)  of  the  batteries  is  86.8%.  The 
minimum  SOC  to  start  the  PEMFC  is  set  to  85%. 


Fig.  6.  NESLAB  RTE  thermostatic  bath  circulator. 

As  the  experiment  is  running,  the  EMS  must  estimate  in  real¬ 
time  the  time  ts t  for  the  SOC  to  reach  85%,  using  only  the  past 
data  from  Pm  and  Eq.  (7). 

7.3.3.  PEMFC  thermal  controller  parameters 

The  fuel  cell  stack  has  no  thermal  insulation  and  the  ambient 
temperature  measured  is  297  K.  To  test  the  different  stages 
involved  in  the  fuel  cell  cold  start  problem  (warming  stage,  self¬ 
heating  stage  and  high  temperature  stage),  we  select  the 
following  operating  conditions  which  can  be  easily  managed  by  our 
test  bench  system: 

•  At  the  beginning  (t  =  0),  the  temperature  of  the  stack  is  the 
same  as  the  ambient  temperature  (Tamb  =  297.15  I<). 

•  The  PEMFC  operating  range  is  set  to  [303.15  K,  323.15  K].  The 
minimum  operating  temperature  Tf™in  =  303.15  I<  is  greater 
than  the  ambient  temperature  in  order  to  allow  the  warm-up 
system  (the  thermostatic  bath)  to  heat  the  stack  and  increase 
TfC  from  Tf°  to  Tf™in.  This  process  will  emulate  the  warming-up 
stage  of  the  cold  start  strategy. 

•  The  maximum  power  for  the  warm-up  system  is  700  W.  The 
cooling  fluid  is  the  water  circulating  in  a  closed-loop  circuit. 

•  From  303.15  K  to  T™ax  =  323.15  K,  the  fuel  cell  is  in  the  self¬ 
heating  phase  where  the  reaction  heat  is  used  to  increase 
quickly  the  temperature  of  the  stack. 

•  As  the  temperature  reaches  Tf™ax,  the  process  regulation  is 
started. 

7.4.  Experimental  and  simulation  results 

In  this  section,  we  present  a  comparative  study  between  the 
simulation  and  the  experimental  results  when  the  test  scenario 
described  above  is  executed  on  our  test  bench. 
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Fig.  7.  User  interface  for  the  monitoring  and  the  control  of  the  PEMFC  system. 


7 A  A.  Validation  of  the  time  to  start  the  PEMFC 

Using  the  power  demand  represented  in  Fig.  8,  the  time  to  start 
the  fuel  cell  is  estimated  based  on  the  history  of  Pm.  Fig.  9  shows  the 
real-time  estimation  of  ts t  which  is  obtained  by  using  Eq.  (7).  For 
comparison  purpose,  we  determined  off-line  the  time  for  the  SOC 
to  reach  85%  by  considering  all  the  time  sequence  of  the  power 
demand.  This  evaluation  (tst  =  1000  s)  is  represented  by  the  dash 
line  in  Fig.  9. 

One  can  observe  in  Fig.  9  that  at  the  beginning,  the  estimation  tst 
was  too  far  from  the  value  of  tst.  As  the  experiment  is  running,  and 
the  number  of  past  data  of  Pm  is  increasing,  the  estimation  error 
becomes  small,  despite  the  big  variation  of  the  power  demand 
signal.  In  Fig.  10,  the  fuel  cell  starts  the  power  production  at 
t  =  1029  s  which  is  greater  than  tst.  This  small  lag  is  due  to:  (i)  the 
stack  temperature  sensor  position  which  measures  a  local 


temperature  within  the  stack  and  (ii)  the  valves  activation  delays 
(anode  and  cathode  valves).  Nevertheless,  when  the  fuel  cell  was 
started,  the  SOC  decreasing  rate  is  significantly  reduced.  This  result 
suggests  that  the  proposed  method  can  be  used  to  perform  a  fairly 
good  real-time  estimation  of  the  time  to  start  running  the  PEMFC. 

7.4.2.  Estimation  of  the  time  to  start  heating  the  stack 

Given  tst,  the  fuel  cell  temperature  controller  must  determine 
the  time  to  start  the  warm-up  system  tst,he  based  on  the  different 
heat  transfer  characteristics,  the  ambient  temperature  Tam b,  the  fuel 
cell  minimum  T£?in,  the  heating  duration  the,  the  time  to  start  the 
fuel  cell  tst  and  the  heating  power  that  can  be  provided  by  the 


Fig.  8.  FC-PHEV  power  demand  based  on  a  down  scaled  experimental  driving  cycle. 


Fig.  9.  Real-time  estimation  of  tst. 
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Fig.  10.  Energy  flow  during  the  experiment. 
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Fig.  11.  Experimental  result. 


warm-up  system  Phe.  As  we  said  before,  lfc  =  Tamb  at  the  beginning 
of  the  experiments. 

Using  Eq.  (18),  the  time  to  start  heating  is  estimated  to  628  s  as 
shown  in  the  lower  graph  of  Fig.  11.  If  this  estimation  is  good,  the 
stack  temperature  should  be  close  to  the  minimum  temperature 
rf™in  at  t  =  tst  so  that  the  fuel  cell  can  start  on  time. 

In  Fig.  11  (first  graph),  one  can  observe  that  the  fuel  cell 
temperature  is  close  to  30  °C  which  corresponds  to 
Tf™in  =  303.15  K  at  t  =  1029  s.  Although  the  presence  of  this  lag  of 
29  s  (which  represents  only  2.9%  of  tst  value),  this  experimental 
result  suggests  that  the  method  to  estimate  the  time  to  start 
heating  the  fuel  cell  is  good  enough. 

7.4.3.  Fuel  cell  thermal  model  during  the  stack  self  heating  stage 

As  the  fuel  cell  temperature  reaches  Tf™in,  the  fuel  cell  is  auto¬ 
matically  started  (see  the  pink  curve  of  the  second  graph  in  Fig.  11 ). 
The  self  produced  heat  by  the  stack  is  shown  in  the  second  graph, 
the  green  curve,  whilst  the  electric  power  produced  by  the  fuel  cell 
is  represented  in  the  same  graph,  the  pink  curve.  The  stack 
temperature  rises  from  24  °C  (297.15  I<)  to 
50  °C(Tf™ax  =  323.15 1<).  During  the  self-heating  stage,  the 
comparison  between  the  simulated  fuel  cell  temperature  and  the 
measured  temperature  indicates  that  both  curves  have  a  similar 
trend.  However,  the  experimental  temperature  is  shifted  down  by 
7%  because  of: 

•  The  fuel  cell  nonlinear  temperature  behavior  that  is  not  well 
included  into  the  developed  lumped  model  (see  Eq.  (22)). 


•  The  position  of  the  temperature  sensor  in  the  stack.  A  more 
interesting  way  to  capture  the  fuel  cell  average  temperature 
would  be  to  use  several  sensors  at  several  positions  within  the 
stack. 

For  automotive  application  with  a  serial  topology,  this  level  of 
mismatch  between  the  lumped  model  and  the  experimental  one  is 
acceptable  since  there  is  no  need  for  the  fuel  cell  to  precisely  track 
the  power  demand. 

Since  the  maximum  heating  power  was  used  by  the  warm-up 
system  (see  Fig.  11,  the  second  graph,  the  blue  curve),  the  stack 
heating  duration  is  the  shortest  one.  Otherwise,  a  more  powerful 
heating  system  should  be  used. 

7.4.4.  Global  strategy  simulation 

The  simulation  aims  at:  (i)  providing  the  stack  thermal  behavior 
when  the  starting  temperature  belongs  to  the  sub-freezing  range; 
(ii)  comparing  the  energy  used  to  warm-up  the  stack  during  a  given 
trip.  We  consider  an  extended  version  of  the  vehicle  power  demand 
used  for  the  experimentation  and  shown  in  Fig.  8.  The  extension 
power  demand  profile,  obtained  by  concatenating  several  times  the 
experimental  one,  has  a  duration  of  5000  s.  In  addition,  we  suppose 
that: 

•  the  initial  stack  temperature  is  -20  °C  (253.15  I<); 

•  the  stack  maximum  operating  temperature  is  T£}ax  =  60  °C 
(233.15  I<); 

•  the  minimum  operating  temperature  is  T£}m  =  5  °C. 
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Fig.  12.  Simulations  results. 


During  the  warm-up  phase  (when  the  stack  temperature 
belongs  to  [-20  °C,  5  °C]),  we  used  the  same  external  warm-up 
system  characteristics  as  the  experimentation  one.  However,  in 
order  to  show  the  heat  losses  through  the  natural  convection,  we 
analyze  two  scenarios.  In  the  first  one  (called  scenario  1),  we 
suppose  that  the  selected  warm-up  system  operates  at  its 
maximum  heating  power  (700  W).  This  scenario  corresponds  to  the 
heating  strategy  proposed  in  this  paper.  In  the  second  scenario 
(called  scenario  2),  we  consider  that  the  warm-up  system  starts 
heating  with  a  constant  heating  power  which  is  less  than  the 
maximum  one  and  which  will  let  the  stack  temperature  to  reach 
Tf™m  at  the  prescribed  time  to  start  the  fuel  cell  (fst).  During  the  self¬ 
heating  phase,  the  proposed  time  optimal  strategy  is  used. 

Fig.  12  shows  the  simulation  result.  The  graphs  (a)  and  (b)  are 
related  to  scenario  1  whereas  the  graphs  (c)  and  (d)  are  related  to 
scenario  2.  With  both  scenarios,  the  fuel  cell  temperature  reaches 
the  Tf™in  at  the  prescribed  time  ts t.  However,  during  the  warm-up 
phase,  the  heat  losses  due  to  natural  convection  with  scenario  1 
(green  curve  in  the  graph  (a))  is  less  than  convection  losses  with 
scenario  2  (green  curve  in  the  graph  (c)).  The  energy  used  by  the 
warm-up  system  with  scenario  1  and  with  scenario  2  are  0.304  MJ 
and  0.317  MJ,  respectively.  These  results  indicate  that  the  proposed 
heating  strategy  in  the  warm-up  phase  is  effective. 

8.  Conclusion 

A  global  strategy  for  a  PEMFC  cold  startup  is  investigated 
through  a  theoretical  analysis  and  an  experimental  validation  on 
a  laboratory  test  bench.  This  strategy  aims  at  providing  an  efficient 
method  to  manage  a  PEMFC  thermal  behavior  within  the  frame¬ 
work  of  a  Fuel  Cell  Plug-in  Hybrid  Vehicle  (FC-PHEV)  with  a  serial 
topology.  The  Energy  Management  System  (EMS)  based  on 
a  supervisory  architecture  is  proposed,  implemented  and  validated. 
This  system  estimates  in  advance  and  in  real-time,  the  time  to  start 
the  fuel  cell  based  on  the  time  for  the  battery  energy  to  reach 
a  minimum  prescribed  energy.  Given  this  estimation  and  the 
amount  of  additional  energy  required,  the  fuel  cell  temperature 


management  strategy  computes  the  most  appropriate  time  to  start 
heating  the  stack  in  order  to  reduce  the  heat  loss  through  the 
natural  convection.  To  quickly  increase  the  temperature  when  the 
fuel  cell  is  started,  we  propose  the  use  of  the  reaction  heat  as  a  heat 
power  source  through  a  time  optimal  control  law.  The  experimental 
results  have  shown  that  the  proposed  approach  is  efficient  and  can 
be  implemented  in  real-time  on  FC-PHEVs. 
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Nomenclature 

Alphabets 
c:  specific  heat 
C:  specific  heat  capacity 
E:  energy 

F:  Faraday  constant 

AH:  enthalpy  of  formation 

h:  convection  heat  transfer  coefficient 

I:  electric  current 

m:  thermal  mass 

m:  mass  flow  rate 

P:  power 

Q:  rate  of  the  heat 

S:  effective  heat  transfer  surface 

t:  time 


T:  temperature 
V:  voltage 

Greek  symbols 
ry.  efficiency 

Xa:  stoichiometric  air  ratio 

Superscripts  and  subscripts 
a:  air 

act:  activation 

amb:  ambient 

b:  battery  pack 

cell:  cell  of  the  stack 

con:  concentration 

cw:  water  forced  convection 

fc:  fuel  cell 

he:  warm-up  system  heater  element 

in:  inlet 

m:  mechanic 

max:  maximum 

min:  minimum 

mst:  mass  species  transportation 
nc:  natural  convection 
ohm:  Ohmic 
out:  outlet 

ocv:  open  circuit  voltage 

rh:  hydrogen  with  oxygen  reaction 

st:  start  time 

td:  road  trip  duration 

tr:  traction 

u:  vehicle  user 

w:  fluid 

whe:  fluid  heat  exchanger 


